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We discuss the issue of interplay between perturbative and non-perturbative phenomena for power corrections 
to e + e~ event shapes. 
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1. Introduction 

For some time now it has been known that 
to describe event-shape measures in e + e _ colli- 
sions, perturbative QCD on its own is not suf- 
ficient. In general there is a need for signifi- 
cant additional corrections, phenomenologically 
of the same size as the next-to-leading correc- 
tions, of non-perturbative (NP) origin. Initial 
estimates for these corrections came from Monte 
Carlo event generators |]J| . More recently, inter- 
est has arisen in analytic approaches , based 
on ideas such as the high-order behaviour of per- 
turbation theory (renormalons), or the concept of 
an infra-red-finite coupling, which predict correc- 
tions of order l/Q, where Q is the centre-of-mass 
energy. This form is in rather good agreement 
with the experimentally required corrections, but 
there is no a priori method (except perhaps from 
the lattice) of determining the coefficient in front 
of l/Q: it is a fundamentally non-perturbative 
quantity. There is however the possibility that 
relative coefficients, from one observable to the 
next, may be estimated. This is referred to as 
the universality hypothesis. 

As will be discussed in section |^, traditional 
methods for calculating the relative coefficients 
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involve considering a qq pair together with a sin- 
gle soft gluon. It turns out that event shapes can 
be divided into two categories according to their 
behaviour in the presence of a qq pair and many 
soft gluons: those whose value is a linear com- 
bination of contributions from each of the soft 
gluons, and those which instead involve a non- 
linear combination. The traditional methods are 
adequate for linear observables. But they are not 
suitable for non-linear ones because the presence 
of perturbative (PT) soft gluons requires that one 
take into account correlations between PT and 
NP gluons. Section ^ illustrates how this is done. 

1.1. Event-shape definitions 

Let us first define the event-shapes that will be 
considered here. The thrust is given by, 



T 



(i) 



and measures the extent to which an event is 
'pencil-like'. The C-parameter is 
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The heavy jet mass is, 
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where M% is the invariant mass of the heavier of 
the two hemispheres defined by a plane perpen- 
dicular to the thrust axis n. Finally we have the 



total jet-broadening, 



(4) 



where the k t are measured perpendicular to the 
thrust axis. A similar quantity, known as the 
wide-jet broadening, Byy, is defined as the larger 
of the broadenings calculated separately in the 
two hemispheres. 

These events shapes have the property that in 
the limit of a two-jet event, 1 — T, C, ph, Bt and 
Bw go to zero. 

2. Traditional power corrections 

There are many related approaches to the cal- 
culation of power corrections to event shapes (and 
to a whole range of other observables): renor- 
malons ||, the massive gluon approach J|, the 
dispersive approach ||, the two-loop approach 
Q. In event-shape studies, they essentially re- 
duce to considering the behaviour of the observ- 
able in the presence of the qq pair and a single 
soft gluon. For a soft gluon with transverse mo- 
mentum kt and rapidity 77, the value of the event 
shape variable V is given by kt/Q times a 'char- 
acteristic function' fv(v) : 
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plus terms of order k\jQ 2 . Schematically, the 
power correction can be seen as coming from the 
integral of this characteristic function over rapid- 
ity and transverse momentum, including a factor 

Of Os(fct). 

More specifically, one considers only the non- 
perturbative part of a s , (which we'll call cknp), 
which is non-zero only for small scales of the order 
of a GeV. So we have an expression for the power 
correction which isQ 

2C F f dk t k t 



SV = 



k t Q 



"Np(^t) / dry f v (77). 
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2 In practice one includes an additional factor of 2M/n ~ 
1.14, related to the splitting of the soft-gluon and details 



Since all the observables have the same depen- 
dence on kt, the integral over kt yields the same 
value in each case: 



dk t kt . . 



w „ a s {kt) - opT(fct) 
dk t — 



Q 



(6) 

where apT is the perturbative expansion for a s 
around scale Q; the integral can be truncated 
at [ij ~ 2 GeV, because beyond that scale cjnp 
should die off very quickly. 

The observable-dependent part of the power 
correction comes from the integral over rapidity, 
which yields an observable-dependent coefficient, 
cv, 



c-v = I dry f v (77) , 



(7) 



multiplying a^/Q. 

The postulate that ao experimentally really is 
the same for all variables, is referred to as univer- 
sality. A priori we have no way of calculating ao 
(except perhaps on the lattice ||), but we should 
be able to measure it in different observables and 
find it to be consistently the same. 

3. Interplay with the perturbative event 

The power correction as calculated above is 
valid for a situation in which the non-perturbative 
gluon which we consider (termed 'gluer' by Dok- 
shitzer ) is the only particle present apart from 
the initial qq pair. In practice however the event 
nearly always contains many soft and collinear 
gluons of perturbative origin. There could also 
be several 'gluers'. 

So we need to understand what happens to the 
event shape in the presence of more than one soft 
particle. 

3.1. The thrust and C-parameter 

Some variables, such as the thrust are simple. 
In the presence of many soft particles, 1 — T just 
becomes 



of the definition of a s (kt) 



1-T: 



-toil 
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Figure 1. Situations of relevance for the heavy-jet 
mass. Solid straight lines represent the underly- 
ing perturbative event. Curly lines are gluers. 
Dashed gluers do not contribute to the heavy-jet 
mass. 



i.e. it is a linear combination of the contributions 
from each soft particle. Thus the contribution of 
the gluer will just add to that of the underlying 
perturbative event, and the consideration of the 
power correction based on the presence of just a 
single gluon will remain valid. A similar argu- 
ment applies to the C-parameter. 

3.2. The heavy-jet mass 

Considering the heavy-jet mass ph, we find that 
the situation is somewhat different. As given 
above, in the presence of a single gluon (figs. |l|a 
and |l|b), ph is just k t /Qe~^, because the heavy- 
jet is always the one containing the gluon. In the 
presence of two particles in the same hemisphere 
(fig. |l|c), then the contribution is 

I(fe t)ie -M+fc ti2 e-tel) , (9) 

which is fine because it is a linear combination of 
the two contributions. But when the two glu- 
ons are in opposite hemispheres (fig. |l|d), the 
heavy jet mass is the larger of the two hemisphere 



masses: 

imax^e-^U^e- 1 " 21 ) , (10) 

This is a non-linear combination of the two con- 
tributions, which means that in the presence of 
the full PT event, the power correction cannot 
simply be estimated by adding the 1-gluon result 
to the perturbative answer. 

So how do we determine the power correction 
in this case? The fundamental point is to ob- 
serve that the contribution to the event shape 
from the underlying perturbative event is always 
much larger than that from the gluers. To see 
why, consider that because of Sudakov suppres- 
sion of events without soft and collinear gluons, 
the probability of an event-shape variable V be- 
ing smaller than some value V max <C 1 is roughly 

P(V < V max ) ~ exp [~N^- In 2 V ma ^j 

(11) 

where N is a variable-dependent number. Thus 
typical values for V are of the order of 



^typical ~ exp ( - -— J - exp ( - y/lnQ 



The contribution from non-perturbative effects, 

Vnp ~ ~ exp (- InQ) , 

is much smaller. 

As a result, in the limit of large Q, it will be the 
perturbative event which 'decides' which hemi- 
sphere is heavy (figs. [l|s and [j]f). Once that infor- 
mation is known, ph is just a linear combination 
of the contributions from all the particles in the 
heavier hemisphere. So by invoking the large dif- 
ference between the typical perturbative and non- 
perturbative scales, we have been able to reduce 
the heavy-jet mass to being linear in the contri- 
butions from the gluers. Since it is only gluers in 
one of the hemispheres that are relevant, we halve 
the power correction compared to the single gluon 
estimate p0|j7[. 
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Figure 2. Broadening for different quark axes. 



3.3. The jet broadenings 

The jet-broadenings give a particularly rich 
example of a non-linear event-shape observable 
which is reduced to being linear in the presence 
of perturbative radiation. 

In the presence of a single soft gluon there is 
a contribution from the transverse momentum of 
the gluon and from the recoil momentum of the 
quark. Noting the factor of a half in the defini- 
tion of the broadening, eq. (0), one obtains that 
the contribution to B is kt/Q- In the presence of 
two soft gluons in the same hemisphere, the sit- 
uation becomes more complex: the contributions 
to B from the k t 's of the soft gluons will simply 
be (kt i + fc t ,2)/2Q; but the contribution from the 
quark recoil will be |fct t i + fet,2|/2Q, which is any- 
thing but simple, especially if one starts having to 
take into account many gluons. This first compli- 
cation is resolved by noting that recoil from per- 
turbative radiation will cause the quark to have 
a large (relative to the NP scale) transverse mo- 
mentum p t . After carrying out the integration 
over the azimuthal angle of the gluer one obtains 
that the recoil due to the gluer is 

J ^ (Vp? + *? + 2k t p t cos 4>-p\~h- 

(12) 

Being quadratic in kt, this will lead at most to 
a 1/Q 2 correction. So in the presence of PT ra- 
diation we can neglect quark recoil. This halves 
one's expectation for the power correction (since 
in the 1-gluon approximation equal contributions 
came the gluon and from the quark recoil), but 



still leaves it proportional to (lnQ/A)/Q. This 
form arises because the integral over rapidity, (0) 
is bounded by the kinematic limit which is ap- 
proximately InQ/fct — InQ/A. Experimentally, 
this was found by both the HI and JADE collab- 
orations to be incompatible with the data 0,^2) . 
The JADE collaboration made the additional ob- 
servation that the power correction seemed to de- 
pend on the value of B itself. 

The further element that has been neglected is 
that of the axis with respect to which one mea- 
sures the transverse momenta. One assumes a 
dkt/kt distribution for emission, where k t is given 
with respect to the quark axis, but the broad- 
ening measures transverse momenta with respect 
to the thrust axis. The extent to which these 
two axes coincide affects the power correction. 
This is illustrated in figure 0. If the quark axis 
and thrust axis coincide, fig. Ha, then gluons with 
(the same k t and) large and small angles to the 
quark axis contribute equally. If on the other 
hand the quark axis does not coincide with the 
thrust axis, fig. ||b, then gluons at large angles to 
the quark axis contribute as before. But gluons 
close to the quark axis (angles smaller than the 
angle between the thrust and quark axes) do not 
contribute — or more precisely their contribution 
is exactly cancelled by a corresponding longitudi- 
nal recoil of the quark. Therefore the rapidity 
integral for the gluer does not extend beyond the 
quark rapidity, r) q , leading to a power correction 
proportional to rj q /Q. To determine r] q , one ob- 
serves that it is determined by the hardest gluon 
(transverse momentum p t ) in the event, so that 
r\ q ~ hiQ/pt ~ hxl/B. Thus we arrive at the 
result that (for a single hemisphere — Bi), the 
coefficient of the power correction is jll| 

c Bl { Bl )= l -\n^ + 0{l) . (13) 

This applies separately for the broadenings from 
each hemisphere. 

Mean Broadenings 

Let us first use this information to determine 
the power correction to the mean broadenings. 
We have to integrate the power correction with 
the perturbative distribution for the broadening. 
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Figure 3. Fit to the mean jet-broadening, using 
a power correction of the form The lowest 

line is the LO prediction, the next line includes 
the NLO corrections, and the highest line includes 
additionally the power correction. 



This has to be done separately for each hemi- 
sphere. The perturbative distribution for the 
single-hemisphere broadening goes roughly as 



B x da 2a s C F , 1 

~ In — exp 

cr dBi 7T Si 
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(14) 



Integrating this with ( |l3| ) and including a factor 
of two to take into account the two hemispheres 
gives us the coefficient of the power correction for 
the broadening: 



■0(1), 



(15) 



where the terms of (1) are to be found in J13[ ]. 
This form for the power correction is used in fig- 
ure H where one sees good agreement with the 
data. 



For the wide-jet broadening, the distribution is 
similar to ([l4| ) but with an extra factor of two in 
front of each a s . This together with the the fact 
that we have the NP contribution from just one 
hemisphere, gives a power correction coefficient 
of 



•0(1) 



(16) 



Broadening distributions 

The situation for the distribution of the total 
jet-broadening is a little trickier. Essentially, for 
small jet-broadenings each hemisphere contains 
about half the total broadening and the power 
correction is roughly twice (|l3|): 



c BT (B T )~]n— + 0(1) . 

JDJ 1 



(17) 



For larger jet-broadenings, most of the jet- 
broadening comes from one hemisphere (whose 
power correction is just (]l3|)), while the other 
hemisphere is essentially free to have almost any 
value of the broadening and so its power correc- 
tion is half of (R5h giving a total of 



c Bt {B t ) 



2\/o. s Cf 



0(1) 



(18) 



Again, the full formulae, which interpolate be- 
tween the two regimes, are given in |l3| ]. Two ex- 
amples of the broadening distribution compared 
to data are shown in figure || and the need for and 
success of the power correction are clearly visible. 

The power correction to the distribution of the 
wide-jet broadening is much simpler since we only 
have one hemisphere to worry about. It is given 
simply by (|l5|). The comparison to the data is 
however quite unsatisfactory. The distribution is 
shown together with OPAL data []14| in figure [|. 
There is quite clear disagreement, with a substan- 
tial need for the whole distribution to be squeezed 
to smaller B\y, especially at moderate and large 
values of Bw- Currently the origin of this prob- 
lem is not understood. But it should be noted 
that the wide-jet broadening is in some respects 
quite a subtle variable: for example, going from a 
three-jet to a four-jet event can reduce the value 




Figure 4. The distribution of the total jet-broadening from the JADE |l2| and OPAL jl4| collaborations, 
compared to theoretical predictions: the dashed lines are the resummed perturbative prediction, while 
the solid lines include the power correction. The parameters, a s = 0.1158 and o?o = 0.5368, have been 
obtained from a fit to Bt distributions at a range of energies. 



of Bw- These subtleties may translate into large 
higher-order corrections, both to the perturbative 
and non-perturbative parts. 



4. Conclusions 

Figure ^ shows 2-a contours from fits to a va- 
riety of observables, both for distributions and 
mean values. Given that one does not necessar- 
ily expect ceo to be consistent between one ob- 
servable and another to better than 20%, because 
of potential higher-order non-perturbative effects, 
and that theoretical systematic errors have not 
been taken into account, the agreement for a s 
and «o between observables is fairly reasonable. 
The only exception is the wide-jet broadening, 
which as we have seen above has yet to be fully 
understood. 

Thus the approach of power corrections is a 
powerful tool in the analysis of event shapes. An 
important element in its success is a treatment 
of the interplay between perturbative and non- 
perturbative radiation. 
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